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Abstract Several path integral representations for the 
T- matrix in nonrelativistic potential scattering are given 
which produce the complete Born series when expanded 
to all orders and the eikonal approximation if the quan- 
tum fluctuations are suppressed. They are obtained with 
the help of "phantom" degrees of freedom which take 
away explicit phases that diverge for asymptotic times. 
Energy conservation is enforced by imposing a Faddeev- 
Popov-like constraint in the velocity path integral. An 
attempt is made to evaluate stochastically the real-time 
path integral for potential scattering and generaliza- 
tions to relativistic scattering are discussed. 

Keywords Path integrals • Potential scattering • 
Real-time Monte Carlo method 
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1 Introduction 

Nonrelativistic quantum mechanical scattering in a lo- 
cal potential is usually described in the framework of 
time-dependent or time-independent solutions of the 
Schrodingcr equation. Path-integral methods in quan- 
tum mechanics, on the other hand, are mostly applied 
to the discrete spectrum, e.g. for harmonic or anhar- 
monic oscillators. In contrast, the transition matrix for 
the continuous spectrum is rarely represented as path 
integral. Even if available, many representations turn 
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out to be rather formal, e.g. requiring infinitely many 
differentiations or infinite time limits to be performed. 
This is not only impractible but also unfortunate since 
a convenient path integral representation may lead to 
new approximations and may be extended readily to the 
many-body problem or Quantum Field Theory. Also 
the long-standing problem how to evaluate real-time 
path integrals by stochastic methods needs a suitable 
path integral representation as starting method. There 
has been significant progress in dealing with real-time 
path integrals for dissipative systems but in closed sys- 
tems and infinite scattering times only zero-energy scat- 
tering seems to be tractable by Euclidean Monte-Carlo 
methods at present pQ. 

Here I will describe an approach developed recently 
[5] to overcome these difficulties while still being prac- 
tical. Its main features are the use of "phantom" de- 
grees of freedom to cancel phases which would diverge 
for asymptotic scattering times. The eikonal approx- 
imation - valid for high energy and small scattering 
angles and the basis for Glauber's multiple scattering 
approach - is "built in" from the beginning as it is ob- 
tained when all quantum fluctuations are neglected. 



2 How to obtain a path-integral representation 
of the T-matrix 

We start with the definition of the 5-matrix as matrix 
element of the time-evolution operator in the interac- 
tion picture at infinite scattering times 



S^ f = lim U f U!(T,-T) 

T— too \ 

= lim e^ Ei+E ^ T U f 



U(T, -T) 
2niS(E l -E f ) Ti^f (1) 
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and its connection with the T-matrix. Here the scat- 
tering states are normalized according to (4>f\4>i) — 
(27r) 3 (5 (3) (k / - k,) . With the definition of Vi in terms 
of the full time-evolution operator U one could use its 
standard representation as a path integral over trajec- 
tories (ft, = 1 ) 

,.x(T)=x b 

U(x b ,T;x a ,-T) = / V 3 x(t) 

Jx(-T)=x Q 



x exp 



1 J**&(j*(t)*-V(x(t)) 



(2) 



connecting initial and final points. However, keeping the 
boundary conditions of the paths is cumbersome and it 
is more convenient to convert the path integral into one 
over velocities by inserting 



1 = / V 3 v5 {3) (v(t) -x(t)) . 
Then one obtains 
(5-1)-^, = lim e**^ [d 3 re- 

T^oo J 

r+T 



(3) 



V v exp 



1 ! T dt j v W 




/" +T / K 
i / cftV" r H t + x (t) 

7-T V m 



- 1 



,(4) 



x v (t) 
Here 



+T 



dt'sgn(t-t') v(i'), x^(t) = u(i).(5) 



q = k/ ki , K = 



ki + k/ 



and the path integral is normalized such that 



V v exp 



1 ! dt j v w 



i 



(6) 



(7) 



There remain two problems to obtain the desired path- 
integral representation of the T-matrix: 



1. There is still an explicit phase 
<2>(T) 



K 2 \ q 2 



(8) 



which seems to produce a divergence in the limit 
T — > oo. Of course, it is cancelled in each order of 
perturbation theory as one can check. 
2. Energy conservation is not yet evident. 

1 Better done in the discretized version of the path integral 
which shows that one has ./Vv-integrations compared to N—l ones 
over the intermediate points of the trajectory when the endpoints 
are fixed. 



The first problem is solved by realizing that this 
phase factor is generated by applying exp(— iT Aj (4m)) 
on exp(— iq ■ r). Integrating by parts and "undoing the 
square" by means of a path integral over an "antiveloc- 
ity" 

-TA\= I V 3 w exp 



exp 



i 

Am 



—T 



, ™ 9 / x 



± 



+ T 



1 



— T 



dt-f(t)^{t) 



V 



(9) 



this simply becomes a shift operator. In Eq. the arbi- 
trary real function only has to fulfill J dtf 2 (t) = 2T 
and we choose it for simplicity as f{t) = sgn(i). Note 
that the sign of the kinetic energy for the antivelocity 
is opposite to the usual kinetic energy in order to ob- 
tain a real shift in the argument of the potential. This 
has an amusing similarity with the "phantom" degrees 
of freedom in the Lee- Wick approach to quantum field 
theory [3] which has been discussed again recently [3]. 

Without a perturbative expansion of the S'-matrix 
the second problem is solved by the Faddeev-Popov 
trick of multiplying the path integral with 



1 



IKI 



+ 00 



drS K 



K 



(10) 



and shifting the co-ordinates. Arguing that in the limit 
T — > oo the action is invariant under a finite time shift 
r one obtains 

r+oo / K 

-iq • — r 



= 2tt<5 




d A be 



iq-b 



— iqb 



..(ii) 



Actually, this procedure is a delicate exchange of time- 
limits but it has been checked that the outcome is cor- 
rect by deriving the Born approximation to all orders 
(see appendix of Ref. 0). So we obtain the following 
path-integral representation of the T-matrix 



V 



(3-3) 



K 



d 2 6e -iq-b / V 3 vV 3 w 



x exp 



, m 
dt — 
2 



v 2 (i) 



! (*) 



: j exp[ix(b,v(i),w(i)) ] - 1 



(12) 



with 

X(b,v,w) 



dtV(Z(t)) 



K 



€(t) = b+ — t + x v (t) -x w (0) 



(13) 
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Fig. 1 (Color online). Scattering geometry for a potential of ra- 
dius _R showing the impact parameter b, the ray made by the 
incoming and outgoing momenta kj * , and the mean momentum 

K= (k 1 ;+k / )/2. 



Note that = K/m+ v(t) and 
|K| = fccos(6>/2), |q| = 2fcsin(6»/2) 



(14) 



where 9 is the scattering angle and E — k 2 /(2m) the 
scattering energy. Eqs. (IT2"T) and (fT5|) show that the par- 
ticle travels along a straight-line reference trajectory 
parallel to the mean momentum and that the quantum 
fluctuations are taken into account by the functional 
integral over velocity and antivelocity. 

In Eq. (IT21 a 3-dimensional antivelocity is used as 
indicated by the superscript "(3-3)". By a simultane- 
ous shift of velocity variables and impact parameter b 
it is possible to obtain a representation utilizing only 
a 1-dimensional longitudinal (i.e. parallel to the mean 
momentum K) antivelocity: 



(3-1) 



K 



i— l d 2 be- l(l b I V 3 vVw 
m 



x exp < i 



x < exp 



, m 

dt- 



i I dtV(^ y (t)) 



with 

Cray (^) 
Pray(*) 



Pray(i) 



t + x v {t) -x vX (0) -x w (0)K 



K + ^ sgn(t) = k, &{-t) + k f 0(t) . (16) 



Here the particle travels along a ray made by the inci- 
dent momentum for t < and the final momentum for 
t > 0. This is depicted in Fig. [TJ 

What happens for a nonlocal potential V? In this 
case we have to start not from the Lagrangian path 
integral @ but from the more general phase-space path 



integral 

,x(T)=x t 

C/(x 6 ,T;x 0) -T)= / V 3 x(t) / V 3 p(t) 

Jx(-T)=x a 




P (t) • i(t) - 

%(x(i),p(t)) 



P 2 (t) 
2m 



Here 
V w (x,p) 



py 



(17) 



(18) 



is the Wigner transform of the quantum-mechanical 
operator V. Using it in the phase-space path integral 
is equivalent to applying Weyl's quantization rule to 
velocity-dependent interactions (the "mid-point" rule). 
One can now follow all the previous steps to derive 
a path-integral representation of the T-matrix. Writ- 
ing the additional (unconstrained) momentum variable 
p(t) = mu(t) and using a 3-dimensional antivelocity 
we then arrive at the following expression 

=i £ f d 2 be-^ h f V 3 vV 3 wV 3 u 

l_>/ TO J 

x exp < i 



, TO 

dt — 
2 



v^)-w^(t)-u^) 



x | exp [ i x(b, v, w, u) ] — 1 | 
where the phase function(al) is now 



(19) 



X(b,v,w,u) 



dt V w [$(t),K+mv(t)+mu(t)] (20) 



and £(t) has been defined in Eq. p^|) . Again the path 
integrals are normalized such that Gaussian integrals 
are unity as in Eq. ([7]). 



1 > (15) ^ Applications 



3.1 Eikonal expansions 

For high-energy and small scattering angles the parti- 
cle mainly travels along the reference path as shown 
by an appropriate scaling of variables [2]: all quantum 
fluctuations are suppressed by inverse powers of K — 
kcos(9/2) for the representation with a 3-dimensional 
antivelocity and by inverse powers of k for the "ray" 
version. In the former case the leading term is given by 



71 



%K 

TO 



dtV b 



Kt 



-1 ,(21) 
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a version of the eikonal approximation due to Abar- 
banel and Itzykson [5 . Higher-order corrections can be 
calculated systematically and agree with the systematic 
eikonal expansion due to Wallace [6]. In the nonlocal 
case neglecting all quantum fluctuations v(£), w(t), u(t) 
in Eq. (|19[) gives an approximation for the T-matrix 
first derived in Ref. [TJ. 



3.2 Variational approximations 

The new path-integral representations for the T-matrix 
suggest new variational approximations by applying the 
Feynman- Jensen variational principle. This has been 
studied in Ref. [5] with a trial action which is linear 
in the velocities v(i),w(i) and in Ref. [S] with a more 
general (linear + quadratic) trial action. In both cases 
the first correction to the variational result from a cu- 
mulant expansion has also been evaluated and impres- 
sive agreement with partial-wave calculations has been 
found, even for larger scattering angles. 



3.3 Stochastic evaluation of high-energy scattering 

One may try to evaluate the path integral for the T- 
matrix numerically by stochastic methods as is done in 
the Euclidean case. For this purpose (and for simplicity 
in the representation with a 3-dimensional antivelocity) 
we expand the velocity variables in a complete set of 
harmonic oscillator functions 



v(t) 
w(t) 



G 



71=0 



w, 



(22) 



Here to denotes a characteristic time which we take as 
the one needed to traverse the range R of the potential, 
i. e. to = mR/k and the constant C is chosen such that 
the free action becomes 



Writing 



K 



£(t) = b + — i + x quant (i) 
m 

the quantum trajectory is then given by 

t ' 



[In 1 — I V„ 



<7n(0) W„ 



Equant (£) — ^ ^ 
71=0 

where 

Qn(y) ■= [ dx sgn(y - x) u n (x) . 

Z J -co 



(23) 



(24) 



(25) 



Completeness of the harmonic-oscillator wave functions 
implies 



qn [To 



n=0 



<1« 1 - I = — 



T 



\U - to 



(27) 



where the divergence of the infinite sum has been reg- 
ulated by a large scattering time T as before. 

The path integral (fT2")) can now be written as an 
infinite-dimensional integral over the expansion coeffi- 
cients 



77 



(3-3) . 



iK 



d 2 b 



n 

v ti=0 



d 3 v n d 3 w n 



JS 



f| ( fd 3 v n d 3 w n ) e lS « (e* x - l) . (28) 



71=0 

In any numerical calculation the infinite sum over the 
modes n has to be reduced to a finite one involving 
only N modes. Therefore we split the free action into 
two parts , one (Sq) involving the lower modes n — 
0, . . . N — 1 and the other {Sq) depending only on the 
"upper" expansion coefficients n = N, . . . oo. 

For the upper modes we employ the method of "par- 
tial averaging" [TU], i. e. 



oo _ Un (J d 3 v n d 3 w n ) expjiS^ ) e ix 
)n ^ T\n U d3 vnd 3 w n ) exp(iS>) 
~ exp(i(x)w ) ■ 



(29) 



In Fourier path integrals partial averaging has been 
successfully applied to determine the ground-state and 
thermal properties of bound systems (see, e.g. Ref. [TT]V 
It includes the upper modes (up to infinity) at least 
approximately by averaging them with the free action 
thereby improving convergence with the number N of 
explicit modes |12) . 

Here we apply this method to a scattering prob- 
lem. Note that (x)at depends on the lower coefficients 
and acts as a (complex) interaction phase for the 6N- 
dimensional integral over the coefficients v„,w„, n — 
0, . . . N — 1. The free average of the potential over the 
upper coefficients is readily performed as only Gaussian 
integrals have to be evaluated and one obtains 

r + oo 



(X) N = - 



dtV n 



where 



quant 



(*) 



N-l 

E 

71 = 



q„ 



(t) V x quant 



(*)) 



(30) 



(31) 



is the quantum trajectory involving the lower coeffi- 
cients and 



(26) V a{t) (p) ~ V(p) exp 



(32) 
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<7N(t/t ) 



E 



3n(0) 



(33) 



the Gaussian transform of the potential in momentum 
space. The width turns out to be purely imaginary 

t 

M 

where the first term in the bracket comes from the ve- 
locity integration and the last one from the integration 
over the antivelocity. Due to Eq. (|27[) we can write it 

as 

t 



.16 



<TN(t/t ) = — [T- 
Zm 



. N-l 

n=0 



9n(0) 



q " [ T 



(34) 



showing how the arising divergence is exactly cancelled 
by the contribution of the antivelocity. 

Still it is impossible to evaluate the remaining 6N- 
dimensional integral over the lower coefficients without 
any damping. This is provided by giving the particle a 
complex mass H 

mp ar ticic — ► m(l + ir) (35) 
and the phantom the mass m* 

m-phantom > m i 1 ~ iF ) ■ ( 36 ) 

If this is done in the full path integral, r — > oo would 
lead to the eikonal approximation because all quantum 
fluctuations are suppressed in that limit and r — > 
would give the exact result. This supports the hope 
that numerical evaluations at finite r and an extrapo- 
lation of these results to r — > would yield reasonably 
accurate results for the T-matrix. 

However, it is not possible to modify the path inte- 
gral over the upper coefficients in this way as the phan- 
tom must have the same mass as the particle: If it would 
have a different mass, say Am then the second term on 
the r.h.s of Eq. (|3"4")l would be multiplied by 1/A and the 
diverging term T (originating from the particle) would 
not be cancelled. So partial averaging and damping do 
not marry (easily) ... 

If they remain separated, i. e. if only the explicit 
modes are damped then the stochastic results only slowly 
converge to the exact ones if r is decreased. This is 
shown in Fig. [2] where scattering from a Gaussian po- 
tential 

V(r) = V e~ r2/R2 (37) 
at a fixed scattering angle is considered for N = 1. 
It is seen that the relative deviation of the stochastic 
scattering amplitude indeed decreases if the damping is 
made smaller but that one also needs more and more 
Monte-Carlo calls to obtain a statistically valid result. 

Although encouraging these results indicate that fur- 
ther work is needed to obtain a practical Monte-Carlo 
scheme for high-energy potential scattering. For a dif- 
ferent approach see Ref. [Ti] , 

2 Actually, this is the method chosen in Ref. I13| for rigorously 
defining real-time Feynman path integrals. 
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Fig. 2 (Color online). The relative deviation \Af/f\ of the 
stochastic scattering amplitude from the exact partial wave re- 
sult at a fixed scattering angle 9 = 40° for scattering from the 
Gaussian potential | [37l) with 2mVoR 2 = —4 at kR = 4. 



4 Relativistic extensions 

To stay closer to the spirit and the aims of this work- 
shop I will now consider scattering under relativistic 
conditions, i.e. truly high-energy scattering. In my opin- 
ion quantum field theory provides the most general, 
consistent - albeit not the most practical - framework to 
describe these phenomena and therefore I will consider 
a scalar theory of particles (<?) with mass M interacting 
by exchange of scalar particles (<ft) of mass m given by 
an interaction Lagrangian 



Ci = 



g\&\4>- (38) 

This generalized Wick-Cutkosky model [T5HTB] mim- 
ics meson-exchange between massive nucleons and - 55 
years after its introduction - is still the favorite model H 
of many workers in the field mostly for the relativistic 
bound-state problem (see, e.g. Ref. [H]). 

Let us evaluate the 4-point function of the theory: 
neglecting nucleon loops ("quenched approximation") 
and using the Schwinger representation for the nucleon 
propagator in the presence of the meson field 



[-d 2 - m- 

iMJo 



' +g<f>(x 
dT exp 



i0]" 



-(" 
M v 



d 2 - M 2 +g<j)(x)) T 



(39) 



3 Unfortunately, it also has an unstable ground state which is 
ignored by many authors, even 50 years after this unpleasant fact 
has been proven 1 171 ■ For the present discussion I will join this 
group ... 
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the 4-point function for scattering of two nucleons can 
be written in "worldline" form 

i'OO 

Gi(xi,X2, x 3l x^) = / dTi dT 2 exp [—iM (Ti + T 2 )] 
Jo 



V A y 2 



Vi(+Ti)=x 2 r V2(+T 2 )=x 4 

-D A yi / 

, Vi{-T 1 )=x 1 Jy 2 (-T 2 )=x 3 

x exp ^ iS [yi , y 2 ] + iS int [yi , y 2 ] ^ + ( x\ o x 3 ) . (40) 
Here the free action is given by 



-Ti 



M 



(41) 



while the interaction part is retarded because the mesons 
have been integrated out: 



g 2 2 f +T * [ +Tj 



Ti 



d 4 p 

(27T) 4 



1 



p2 _ m 2 _|_ jQ 



exp [-ip ■ (yi(n) - yj( r 2))] ■ (42) 



This looks (superficially) very similar to a nonrelativis- 
tic description: if we neglect the (i = j)-terms, i. e. the 
self-energy of the nucleons, then we have a 4-dimensional 
analogue of scattering from a Yukawa potential. Still 
there are two different proper times Ti , T 2 over which 
one has to integrate finally and a mass renormalization 
is needed to get rid of the divergencies at small proper 
time - of course, all these features are to be expected 
in a relativistic quantum field theory 0. 

To obtain the relativistic T-matrix one has to Fourier 
transform the 4-point function (PfUl) (easily done by con- 
verting to velocity path integrals), extract the energy- 
momcntum-conserving ^-function and remove the outer 
legs, i.e. multiply with the inverse (full) 2-point func- 
tions of the external particles. While the extraction of 
energy-momentum conservation is obvious in this field- 
theoretic framework due to the translation invariance 
of the action, the amputation of the full external nu- 
cleon propagators is more involved. Using the Bloch- 
Nordsieck approximation it is possible to achieve that 
by functional differentiation techniques (see, e.g. Ch. 
10.1 in Ref. [TI5]) but I expect that introducing phan- 
tom degrees of freedom as in the nonrelativistic case will 
do a better job. This is presently under investigation. 

4 The Wick-Cutkosky model belongs to the simpler class of 
super-renormalizable theories and is, of course, a far cry from the 
field theory of strong interactions, QCD. One could paraphrase 
a remark made by A. Einstein in 1921: "As far as the models 
refer to reality, they are not manageable; and as far as they are 
manageable, they do not refer to reality" . 



5 Summary 

Several path integral representations for the T-matrix 
in potential scattering have been discussed, including 
a new one for nonlocal potentials. They serve to derive 
high-energy expansions, variational approximations and 
are a natural starting point for a stochastic evaluation 
of high-energy potential scattering. Many avenues lie 
open for further development, in particular in the rela- 
tivistic domain. 
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Note added after publication: A method for non-perturbative 
amputation of the 4-point function in the Wick-Cutkosky model 
and in scalar QED has been developed by Barbashov and colla- 
borators [20ll21ll22| . See also Refs. I23ll24l . 
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